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Abstract 

We calculate the decay constant ratios f Vc / fjM and f Vb / fr- In the calculation 
we take into account the mock meson structures of the mesons, as well as the 
difference of the wave functions at origin of the vector and pseudoscalar mesons 
studied by Ahmady and Mendel. We find that the different spin structures of the 
mesons much affect the ratios. We incorporate our results in the prediction of the 
branching ratios of B — > Kr\ c . 
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The decay constants /jm and fy of the vector mesons are given experimentally 
from the decay rates to e + e~, but f Vc and f Vb of the pseudoscalar mesons are lack of 
experimental results. So it is necessary to calculate theoretically the ratios f Vc / fj/f 
and f Vb /fy for various applications, for example, for the prediction of B(B + — > 
K + r? c ) from the experimental result B(B + -> K + J/^p) = [1.02 ± 0.14] x 10~ 3 
[p], If. It has been widely assumed that the decay constants of the vector and 
pseudoscalar mesons are almost the same by considering their wave functions at 
origin to be the same and using the Van Royen-Weisskopf formula. However, 
Ahmady and Mendel H calculated the ratio of the wave functions at origin by 
considering the perturbation caused by the hyperfme splitting Hamiltonian, and 
found that the ratio is significantly different from unity. Then f Vc and f Vb become 
very different from fj/$ and fy. In this paper we will further take into account 
the influence of the mock meson structure 11 [5| of the mesons on the calculation 
of the decay constants. We will consider its effects on f Vc , f Vb , and /ju, fy, which 
are originated from the different mock meson spin structures of the pseudoscalar 
and vector mesons. We find that the effects are very important contrary to the 
expectation that they are not important for the mesons composed of two heavy 
quarks. As a result, f Vc and f Vb become very close to fj/^ and fy again, as we will 
show in this paper. 

The decay constants of pseudoscalar and vector mesons are defined by 

< 0| Qi^Q' \M P {K) >= f P K», < 0| Q^Q' \M V (K, e) >= f v m v e fl . (1) 

The Van Royen-Weisskopf formula || for the decay constants is given by fu = 
\Jl2/rriM |^m(0)|, where tum and ^m(O) are the mass and the wave function at 
origin of the meson respectively. The Van Royen-Weisskopf formula is widely used 
for the calculation of the meson decay constants. This formula is obtained in 
the limit that the spinors of the quarks inside meson are approximated to two- 
component Pauli spinors 0, In this paper we are interested in the J/ip and T 
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families. From the Van Royen-Weisskopf formula, we have 



/ U \ 2 = l^ e (0)| 2 



' m„ e |^(0)| 2 



1.040 



l^ e (0)| 2 

*w(o)l 2 ' 



(2) 



The approximation 1^^(0)1 = |\&j/^,(0)| has been usually used, with which we get 
(fvc/fj/4>) 2 = 1-040. However, recently Ahmady and Mendel calculated |^ c (0)| 2 / 1 ^j/^,(0)| 2 
in their interesting work ||] based on the perturbation theory of quantum mechan- 
ics, and obtained the ratio as 1.4 ±0.1. By incorporating this result in @, Ahmady 
and Mendel obtained [|7| 



Ahmady and Mendel also obtained \^ Vb (0)\ 2 /\^ T (0)\ 2 = 1.16 ± 0.06 in their work 
0. Then by using mx = 9.460 GeV @ and m Vb = 9.445 GeV which is calcu- 
lated from (mj — m Vb ) = im c jm^) 2 (rrij/^p — m Vc ) given by the hyperfine splitting 
Hamiltonian, we obtain 



If we had taken |* % (0)| = |*r(0)|, {f Vb /h) 2 would have been given to be 1.002 



account the fact that the hyperfine splitting Hamiltonian makes the wave function 
at origin of the pseudoscalar meson bigger than that of the vector meson, and then 
have obtained the improved results given by (|3|) and (§) which are much larger 
than 1.040 and 1.002 which would have been given if the wave functions at origin 
had been taken to be equal. 

The purpose of this paper is to improve further the results of @ and (f|) by 
taking into account the mock meson structures of the pseudoscalar and vector 
mesons, which are different by their spin structures, as well as the difference of 
the wave functions at origin which has been studied by Ahmady and Mendel. We 
work in the relativistic mock meson model of Godfrey, Isgur, and Capstick || . 
in which the meson state composed of a heavy quark Q' and a heavy antiquark Q 




(3) 




(4) 



from a similar equation to (0). Therefore Ahmady and Mendel have taken into 
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is represented by 



\M P (0) > = I ——^-p= $(p Q , 



(2tt)3/2 J2E Q ,2E Q VN C 



K(PQ'> C ) 6 1(P<5' c ) - a l(PQ'' c ) 6 t(PO' c )] 1° > ( 5 ) 

in the meson rest frame (where pq> = — Pq) in which we work in this paper, where 
the arrow indicates a state with spin up (down) along a fixed axis and c is the colour 
index which is summed. Whereas we wrote the pseudoscalar meson state in (|5D, we 
can also write the vector meson states in the same way with the spin combinations 
for the vector states, which are given by (TT), 1/a/2(TJ. + it) an d (II)- In (§) 
we adopted the normalization of the creation and annihilation operators given 
by {a(p, s), at(p', s')} = (27r) 3 2E 5 ss '5 3 (p — p'), and then the meson state in (|) 
is normalized by < Mp(0)\Mp(0) >= 2mp 5 3 (0), and also in the same way for 
the vector meson states. We take the momentum wave function $(p) in (fjD as a 
Gaussian wave function 

*(p) = * (r) = (^-^ < 6 > 

where \l/(r) is the conjugate wave function in coordinate space. In Fig. 1, we display 
six different inter-quark potentials of the potential models in Refs. |9|]-[T4j], which 
were obtained by fitting the data of the J/ip and T families (mainly their spectra). 
The mean square radii (r 2 ) 1 ^ 2 of the J/%jj and T mesons are about 2.2 GeV -1 
and 1.3 GeV -1 respectively, which are in the confining long-distant linear 
potential range. Therefore, using the Gaussian wave function in is appropriate. 
In particular, since we will calculate the ratios of the decay constants by using (ffT|) , 
the Gaussian wave function is reliable to use in our following calculations. 

Since we are concerned with the matrix elements in the left hand sides of (0) 
with the meson states in (^|), it is convenient to represent the meson states by the 
matrix-valued representations given by 

m Pa p = - < 0| Q' a Qp \M P {0) >, * Vafi = -<0\Q , a Qp\M v (O)>, (7) 



where a, (3 are spinor indeces. With (|7|), the formulas in ([I]) are written as 

Tr( 7° 75 ^ P ) = fpmp, Tr( 7 "* v ) = / y m y ^. (8) 

If both two quarks inside the meson are static, the spinor combinations of u(0)v(0) 

for the pseudoscalar and vector meson states are given respectively by |16|, [T7| 

1 1+7° , 11+7° 
P(0,0) = -- ? =^7 5 , V{0,0,e) = A (9) 

where the polarization vectors of the vector meson are given by e± = (1/ v2) (0, 1, ±z, 0) 
and £3 = (0, 0, 0, 1). By Lorentz boosting the static spinors we obtain \l/p and 
in (0) as 



(2vr) 3 / 2 J2E^2E^ J2m Q ,{m Q , + E Q ,) J2m Q (m Q + E Q ) 



(10) 

where / = P or V, and Sp and SV are respectively P(0, 0) and V(0, 0, e) in @. By 

incorporating (|I0| ) in @, we obtain the following formula for the decay constants 

of pseudoscalar and vector mesons in the relativistic mock meson model: 

, _ 2y/3 r d 3 p , E Q , + m Q , PQ + m Ni/2 , p 2 

11 JmjJ (2k)*/* lPJ V 2E Q , 2E Q > \ + G/ (P Q , + m Q ,)(E Q + m Q ) 

(11) 

where ap = — 1 and ay = +1/3. We note that the above formulas for the decay 
constants have been derived by taking the four-component spinor into considera- 
tion, and it is reduced to the Van Royen-Weisskopf formula in the two-component 
spinor limit which corresponds to taking the p — > limit in the last two factors of 

When the meson and quark masses are given, f Vc and fj/^ can be calculated 
from fli~ll) for a given value of the parameter j3 in @. We obtained them numerically 
as functions of j3 by using mj/^ = 3.097 GeV, m nc = 2.979 GeV @, and m c = 1.78 



GeV ||18 |. We present the results in Fig. 2. We performed the same calculations 



for f Vb and f T with m T = 9.460 GeV §, m Vb = 9.445 GeV, and m b = 5.17 GeV 
n|, and present the results in Fig. 3. 



Using the experimental values T(J/ip — > e + e ) = 5.27 ± 0.37 keV, T(T — > 
e+e") = 1.32 ± 0.03 keV §, and the formula |I| 

Air a 2 4 1 

T(V -> e + e~) = / y cy, where c J/v , = - and c T = -, (12) 

we get 

/ J/v , = 406 ± 14 MeV, / T = 710 ± 8 MeV. (13) 
Then, from ([TIT) we obtain 



[3 J/% b = 0.644 ± 0.022 GeV, fa = 1.360 ± 0.016 GeV. (14) 

In obtaining the above results, we took the quark masses as m c = 1.30 — 1.85 
GeV and rrib = 4.70 — 5.20 GeV, which cover the quark masses in the six potential 
models in Refs. |9[|-||14||. We note that the errors in ([14]) came from the quark 
mass ranges considered, as well as from the experimental errors of the vector 
meson decay constants in (13). From the results of Ahmady and Mendel ||, 



I^ C (0)I7I^(0)| 2 = 1.4 ±0.1 and |^(0)| 2 /|^ T (0)| 2 = 1.16 ± 0.06, we get 

P V c = Pj/* x (1-4 ± 0.1) 1/3 = 0.721 ± 0.042 GeV, (15) 
/3 Vb = p r x (1.16 ± 0.06) 1/3 = 1.428 ± 0.041 GeV, 

since |^/m(0)| 2 = (Pm/V^) 3 from (|). Then by using these values of (3 Vc and (5 m 
in (|TT|), we obtain 

f Vc = 420 ± 52 MeV, f Vb = 705 ± 27 MeV, (16) 

and then the ratios are given by 

f 2 f 2 

(-p-) =1.06 ±0.14, (p) =0.99 ±0.04, (17) 

J J/i> It 

instead of (|3|) and <Q). We obtained the results (|T7|) by taking into account the 
different mock meson spin structures of the pseudoscalar and vector mesons, as 
well as their different values of the wave functions at origin caused by the hyperfine 
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splitting Hamiltonian. We note that the difference between the results in (|TT| ) and 
those in (j^) and (^) has come from the relativistic correction (obtained by taking 
the four-component spinor into consideration) which is calculated by assuming 
a Gaussian wave function for the quarkonium. The obtained ratios in ( |TTD are 
much close to unity compared to the results of (|]) and (^) which were obtained by 
Ahmady and Mendel ||, [?]] by taking into account only the difference of the wave 
functions at origin. 

The errors in the above results (|14|)-([17]) come from the ranges of the values 
of the quark masses m c and mj which we took as m c = 1.30 — 1.85 GeV and 
m;, = 4.70 — 5.20 GeV, as well as the experimental errors in the values of fj/^ 
and fx in fllTf). In order to see the magnitude of the error in the results induced 
by the sensitivity to the quark masses, we present the results which are obtained 
when we use fixed values of m c and in the calculation: When we use m c = 1.78 



GeV and m b = 5.17 GeV 0, (0) becomes (3j^ = 0.638 ± 0.015 GeV, (3 r = 
1.357 ± 0.013 GeV, (0) /3 Vc = 0.714 ± 0.034 GeV, P Vb = 1.426 ± 0.038 GeV, (|I§ 
f Vc = 424 ± 25 MeV, f Vb = 709 ± 20 MeV, and finally (0) becomes {fjfj^f = 
1.09 ±0.07, {fr,Jhf = 1.00 ±0.04. By comparing these values with those in (H) 
and (|l~7D , we can see the sensitivity of the results to the values of the quark masses 
m c and m^. 

The decays of B — > K J/ip are important for the check of the factorization 
hypothesis and the search of the CP violation phenomena in the B meson decays, 
therefore there have been continuous and intensive experimental improvements on 
their measurements. Recently CLEO reported the new results [|[: 

B(B° -> K° J/ijj) = [1.15 ± 0.23(stat) ± 0.17(syst)] x 10~ 3 , 
B(B° -> (892)° J /if)) = [1.36 ± 0.27(stat) ± 0.22(syst)] x 10" 3 , (18) 
B(B + -> K*(892)+ J /ifj) = [1.58 ± 0.47(stat) ± 0.27(syst)] x 10~ 3 , 
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by incorporating the world average 

B(B + -> K + J/ij}) = [1.02 ± 0.14] x 10~ 3 . (19) 

In connection with B — ► K J/ip, the decays oi B K r] c have been intensively 
studied theoretically 0, |2(], plf| , since the decays are expected to be measured in 
near future and by comparing the two decay modes we can get a lot of valuable 
informations on the hadronic structures. Studying the form factors phenomeno- 
logically in detail, Gourdin et al. obtained the following results pCfl : 

T = nBS?j% =T *&' ' 0957 ~ T - 1 ' 259 ' 

r = m^wim = f * x (^) 2 • 0456 £ f * s °- 872 - (20 > 



By incorporating our obtained decay constant ratios ([T7|) in (|20|), from (|18| ) and 
flI5|) we predict 

5(5+ ^ K + r] c ) = T x [(1.11 ± 0.17) x 10~ 3 ] = [0.90 ~ 1.61] x 10~ 3 , 
B(B° ^ K°r] c ) =f x [(1.25 ± 0.32) x 10~ 3 ] = [0.89 ~ 1.98] x 10" 3 , (21) 
B{B + -> iT(892)+ rj c ) = f * x [(1.72 ± 0.60) x 10~ 3 ] = [0.51 ~ 2.03] x 10~ 3 , 
B(B° -> iT(892)V) =T* x [(1.48 ± 0.39) x 10~ 3 ] = [0.50 ~ 1.63] x 10~ 3 . 

In obtaining fl2~l"|) we combined the errors in (|T7|), (0) and ( |I~9| ) by root mean square. 
For the predictions in ([?Tp we used the range in fl2"U| ) of the values of T and T* 



obtained by Gourdin et al. |[20|| . However, T and T* are very dependent (especially 
T*) on the model for the form factors of (B — > X) and (5 — > X*). Ahmady 
and Mendel applied the heavy quark effective theory and obtained T = 1.12 and 
T* = 0.27 0, whose T is inside, but T* is outside of the range in (p0| ). Particularly, 
Deshpande and Trampetic |2^| emphasized that the value of T* is very much 



dependent on the model adopted and pointed out that the measurements of the 
decay B — > K* r\ c will provide a valuable criterion for the model for the hadronic 
form factors. 
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In conclusion, we calculated the decay constant ratios f v J fj/$ and f Vb / fr by 
taking into account the mock meson structures of the mesons, as well as the differ- 
ence of the wave functions at origin of the pseudoscalar and vector mesons. The 
results have been significantly affected by the different mock meson spin structures 
of the mesons. We incorporated the obtained ratios of the decay constants in the 
prediction of the branching ratios of the B — > K rj c decays. 

Acknowledgements 

This work was supported in part by the Basic Science Research Institute Pro- 
gram, Ministry of Education, Project No. BSRI-95-2414, and in part by Non- 
Directed-Research-Fund, Korea Research Foundation 1996. 



9 



References 

[1] F. Abe et al, CLEO collaboration, Phys. Rev. Lett. 76 (1996) 2015 

[2] Particle Data Group, L. Montanet et al., Phys. Rev. D50 (1994) 1173 

[3] M.R. Ahmady and R.R. Mendel, Phys. Rev. D51 (1995) 141 

[4] R. Van Royen and V.F. Weisskopf, Nuovo Cim. 50 (1967) 617; ibid. 51 (1967) 
583(E) 

[5] S. Godfrey and N. Isgur, Phys. Rev. D32 (1985) 189 
[6] J. Kaplan and J.H. Kiihn, Phys. Lett. B78 (1978) 252 

[7] M.R. Ahmady and R.R. Mendel, Z. Phys. C65 (1995) 263; The convention 
of the decay constant of vector meson in this reference is different from ours, 
by the factor my in Eq. ([!]). The formulas in this paper are written with the 
decay constants in our convention given by Eq. ([I]). 

[8] S. Godfrey, Phys. Rev. D33 (1986) 1391; S. Capstick and S. Godfrey, Phys. 
Rev. D41 (1990) 2856 

[9] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane and T.M. Yan, Phys. Rev. 
D17 (1978) 3090; ibid. D21 (1980) 313(E); ibid. D21 (1980) 203 

[10] K. Hagiwara, A.D. Martin and A.W. Peacock, Z. Phys. C33 (1986) 135 

[11] A. Martin, Phys. Lett. B93 (1980) 338 

[12] A.K. Grant, J.L. Rosner and E. Rynes, Phys. Rev. D47 (1993) 1981 
[13] C. Quigg and J.L. Rosner, Phys. Lett. B71 (1977) 153 
[14] J.L. Richardson, Phys. Lett. B82 (1979) 272 



10 



[15] D.S. Hwang and G.-H. Kim, in the Proceedings of the APCTP Inauguration 
Conference, Seoul, Korea, 1996 

[16] A.F. Falk, H. Georgi, B. Grinstein and M.B. Wise, Nucl. Phys. B343 (1990) 
1; A.F. Falk, ibid. B378 (1992) 79 

[17] M. Neubert, Phys. Rept. 245 (1994) 259; Phys. Rev. D46 (1992) 1076 

[18] D.B. Lichtenberg and R. Roncaglia, in the Proceedings of the Workshop on 
Diquarks II, Torino, Italy, World Scientific 1992 

[19] M. Neubert, V. Rieckert, B. Stech and Q.P. Xu, in Heavy Flavours, edited 
by A.J. Buras and M. Lindner, World Scientific 1992 



[20] M. Gourdin, Y.Y. Keum and X.Y. Pham, Phys. Rev. D52 (1995) 1597; [hep- 
| ph/ 950 1360 



[21] P. Colangelo, C.A. Dominguez and N. Paver, Phys. Lett. B352 (1995) 134 
[22] N.G. Deshpande and J. Trampetic, Phys. Lett. B339 (1994) 270 



11 



Figure Captions 



TM. The radial 



Fig. 1. The inter-quark potentials of the potential models in [] 
distance of the horizontal axis is in the unit of GeV" 1 (1 GeV" 1 = 0.197 fm), and 
the potential energy of the vertical axis is in the unit of GeV. 



Fig. 2. f Vc and fj/$ as functions of the parameter (3. 



Fig. 3. f Vb and fr as functions of the parameter (3. 



12 




p(Gev) 




p(Gev) 



